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ABSTRACT 

This work considers the idea of massive black holes being the constituents of the Galactic 
dark matter halo. It constrains the maximum black hole mass to Mbh ~ 5 x 10 4 M Q by 
examining their influence on the population of globular clusters in our Milky Way. In 
the adopted halo model, globular clusters are exposed to constant bombardment of halo 
objects on their orbits through the Galaxy and thus will steadily gain internal energy. 
Depending on the mass of these halo objects and the structural parameters of the 
globular clusters, they can be disrupted on time scales of a few billion years and below. 
These disruption time scales are calculated using a modification of the well known 
(classical) impulsive approximation and compared with direct N-body simulations of 
such encounter events to ensure the method works correctly. For a set of ten prototypical 
globular cluster models and black hole masses ranging from 10 3 to 10 7 M , Monte- 
Carlo-simulations of 10 000 encounter histories over the period of 10 billion years were 
calculated each, at three different galactocentric distances R = 5, 10 and 15kpc. These 
data were compared with the real globular cluster population in our Galaxy and used 
to obtain the above constraint of Mbh 5 x 10 4 M Q . 
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1 INTRODUCTION 

Comparison of dynamical and photometrical mass determi- 
nations in galaxies like our own strongly indicates the ex- 
istence of a non-luminous dark matter component, that ac- 
counts for most of the mass in these galaxies. The constancy 
of Hi-rotation curves even far beyond the optical edge of 
galactic disks furthermore indicates a l/r 2 -density distri- 
bution of this dark component (Rubin et. al. 1982, 1985). 
Thus one question naturally arises: What is the nature of 
this dark component? A whole variety of dark matter can- 
didates have been proposed in the literature, ranging from 
light leptonic particles to baryonic gaseous matter (possibly 
in the form of cold molecular clouds with fractal structure; 
Pfenninger et. al. 1994a,b) to low mass stars such as white 
dwarfs or "Jupiters" and finally black holes with masses 
Mbh — 1O 6 M0 as remnants of a population of very mas- 
sive primordial stars (for a brief overview see e.g. Chap. 18 
in Peebles 1993 and references therein). 

The latter are of theoretical interest, since the Jeans mass 
at high red shift (just after the epoch of recombination) was 
of the order of 10 6 M (Dicke & Peebles 1968) and black 



holes with this mass could very effectively heat an initially 
cold disk in our Galaxy explaining the observed properties 
(Lacey & Ostriker 1985). But there also are many argu- 
ments against massive black holes as main constituents of 
the Galactic dark halo. Hut & Rees (1992) argued that black 
holes with 10 6 Mq in the inner few kiloparsecs would sink 
into the Galactic center within a fraction of the Hubble 
time due to dynamical friction. This process would lead to 
consequences that strongly disagree with observations. Fur- 
thermore massive black holes would on passages through 
the Galactic disk accrete interstellar matter and thus vi- 
olate background light constraints (Carr 1979). And mas- 
sive black holes would finally cause gravitational lensing of 
cosmologically distant sources. Garrett et al. (1994) stud- 
ied this effect on A 18 cm VLBI maps of the gravitational 
lens system 0957 + 561 A, B and conclude that black holes 
with Mbh > 3 x 10 6 ft -1 Mq cannot account for the dark 
matter in the halo of the main lensing galaxy. Rix & Lake 
(1993) again applied the disk heating scenario to dwarf 
galaxies of the Local Group and derived an upper limit of 
Mbh ^ 1O 4 M . On the other hand, Carr, Bond & Arnett 
(1984) derived a strong lower limit of Mbh > 10 2 Mq by 
examining the possible cosmological consequences of popu- 



2 R. Klessen and A. Burkert 



25 [ 



20 - 



g 
3 



15 - 



10 



0.5 



1.0 1.5 2.0 

c = log 10 (r,/r ) 



2.5 



20 



15 



10 



b) 



50 100 

r, in pc 



150 



10 [ 



Figure 1. 



I 



c) 







1,1 1,1 I, 



d) 



|,| |,| [ 



15 



20 



0123456 5 10 

logioPo in M e /pc a a in km/s 

Distribution of a) concentration index c, b) tidal radius rt, c) logarithm of the central star density log 10 po ar >d d) central 



velocity dispersion <tq for 134 globular clusters in the Milky Way (data from Webbink 1985). 



lation III stars. Their black hole remnants must exceed this 
limit, otherwise the enrichment of the interstellar medium 
would be too pronounced. 

This paper facilitates a different approach: it uses the 
fragility of globular clusters in our Galaxy to constrain Mbh • 
Adopting the hypothesis that the dark halo consists entirely 
of massive black holes, globular clusters are subject to con- 
stant bombardment of these objects on their orbits through 
the Milky Way. The observed properties of globular clus- 
ters thus set strong limits on the mass distribution of these 
black holes. To calculate the response of globular clusters to 
encounters with massive black holes we use a modification 
of the impulsive approximation (Spitzer 1958) and test this 
analytical approach by direct N-body simulations of such 
an event. The excellent agreement of both methods shows 
the validity of the assumptions of the analytical approach 
and allows us to study a large region in parameter space 
since this method consumes much less computer time. We 
obtain analytical values for the average internal energy gain 
of globular clusters due to the interaction with halo objects 
in the Milky Way. Since this energy gain is mostly deter- 
mined by only few central encounters, the scatter around its 
mean value is very large. To account for that scatter and to 
be able to make statistically sound predictions, we adopt ten 
model clusters and compute for every one of them 10 000 en- 
counter histories each at three different galactocentric radii 
(R = 5, 10 and 15 kpc) and for a black hole mass range 



between 10 3 and 10 7 Mq, applying Monte Carlo methods. 
If the increase of internal kinetic energy of the system ex- 
ceeds its binding energy, we assume it dissolves. Thus we 
obtain the fraction of globular clusters, which would survive 
for 10 10 years, depending on internal parameters, such as 
concentration and total mass, their galactocentric distance 
and the assumed black hole mass. Comparing that data with 
observations we constrain the maximum black hole mass to 
M BH = 5 x 10 4 M . 

This value agrees with the mass range obtained by other 
authors using a variety of methods. It especially has to be 
compared with the numbers calculated by Moore (1993). 
He used similar arguments than this paper applies, but did 
not take the huge spread in internal energy gain for dif- 
ferent encounter histories into account. He calculated for 
nine different globular clusters in the outer Galactic halo 
one encounter history each. Therefore his estimate is highly 
uncertain. Taking into account different encounter histories 
of individual globular clusters, the present work stands sta- 
tistically on more solid ground than previous work, leading 
to a more precise estimate for the maximal allowed mass of 
compact dark halo objects. 
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Figure 2. Distribution of total mass Mqc an d concentration index c for 113 Galactic globular clusters of the Webbink (1985). The 
shaded circles are the Mqc an d c values of the ten model clusters in our simulations. 



2 GALACTIC GLOBULAR CLUSTERS 

Galactic globular clusters can very well be fitted by a King 
density profile (King 1966, Illingworth & Illingworth 1976). 
King models are lowered nonsingular isothermal spheres, 
which are generated from the distribution function of an 
isothermal sphere by introducing an energy cut-off . This 
results in a radial cut-off at the tidal radius r\ . King models 
form a sequence that can be parametrized by the concentra- 
tion index c = log 10 (r t /ro). The core radius ro is determined 
by the central density po and the central velocity dispersion 
era: 

A compilation of these parameters for 134 galactic globular 
clusters from the Webbink (1985) data set is shown in Fig. |l| 

According to this data we generated ten model clusters with 
different parameter combinations (see Tab. |l|) to simulate 
the influence of massive compact halo objects on globular 
clusters of different sizes and concentrations. Models #0 to 
#3 have average concentration (c = 1.52) and vary in their 
central dispersion from ao = 3.0kms _1 to g.Okms^ 1 . Models 
#4 to #7 are very little concentrated (c = 0.84) and have 
the same variation in ao as above. Due to their lower concen- 
tration their total mass M is smaller by a factor of roughly 
one third. The sequence #5, #8 and #9 ranges from small to 



large concentrations and has fixed ao. Model #1 represents 
the "mean" of the Galactic globular cluster system and will 
be called the standard cluster throughout this paper. For the 
numerical calculations, each of these clusters is represented 
by 25 000 test particles of equal mass. 

Plotting concentration index c against total mass M, Fig. ^ 
depicts that the Galactic population of globular clusters is 
fairly well sampled by the adopted set of model clusters 
(which are represented by the gray shaded circles). However 
this is not true for the central density po, which in reality 
spans a range from 10 to lO 6 M0/pc 3 (see Fig. [j].c), whereas 
the model clusters have po — 5 x 10 3 Mq/pc 3 . This has to 
be taken into account in the later analysis, when comparing 
the numerical results with observations. Figure [| also in- 
dicates an interesting mass-distance relationship. Globular 
clusters at smaller galactocentric radii tend to be more mas- 
sive than more distant ones. Due to equation |l| this holds 
also for po and a: The inner globular clusters have on aver- 
age higher central densities and larger velocity dispersions. 
Closer to the Galactic center, these globular clusters are 
more severely influenced by tidal shocking during passage 
through the Galactic disk (and bulge in the inner few kpc) 
and encounters with massive halo objects. Therefore they 
have to be more massive and sturdy to survive for a Hubble 
time. 
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Table 1. Properties of the ten model clusters used in the simu- 
lations. The parameters listed are central concentration index c, 
total mass M, core radius ro, tidal radius rt, velocity dispersion 
(jq and density po m the central region of the clusters. 
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3 MODIFIED IMPULSIVE APPROXIMATION 

3.1 Impulsive Approximation with Point Particles 

To calculate the energy and momentum transfer between 
two point particles due to solely gravitational interaction, it 
is advisable to study the problem in the reduced co-ordinate 
system, i.e. one separates the relative motion of the two 
particles from the center of mass motion and follows the 
trajectory of the reduced particle as an one-body problem 
(Landau & Lifschitz 1976, Vol. I). The net velocity change of 
the reduced particle perpendicular and parallel to the initial 
direction is: 
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Here M and m are the masses of the two particles and $ is 
the scattering angle for the encounter. The variables b and 
vo indicate impact parameter and initial relative velocity, 
respectively. Avy always points in the opposite direction of 
the initial velocity vo and Avi lies in the scattering plane. 

The change of absolute velocity of particle m is 
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For the second particle M one just has to interchange the 
masses. 

For large impact parameters, b — > oo, the second t erm of 
the expression in brackets dominates and equations (3.a,b) 
reduce to the classical impulsive approximation, which as- 
sumes a very short encounter time such that one particle 
remains fixed in space whereas the other one moves along a 
straight line with constant velocity: 
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Defining (5 = v b/G(M + m), these equations can be written 
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(classical impulsive approximation) (6) 



3.2 Influence of the Encounter Time 

For the (classical) impulsive approximation to be valid, the 
encounter time r has by far to be the shortest relevant time 
scale in the system under study. The particle can then be 
considered as fixed in space and momentum is transfered 
quasi-instantaneously. Thus, if a periodical oscillation of an 
individual star is superimposed on the linear center of mass 
trajectory of the whole system, the orbital period t or b has 
to be much larger than the time scale r of the external per- 
turbation (r <§C t or b). In the adiabatical limit (r S> t or t), 
the orbital energy of the star is not affected. The perturber 
"sees" the star as being smeared out over the whole peri- 
odic trajectory and hence can only interact with its center 
of mass averaged over the whole period; the orbit acts as a 
rigid body. 
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Figure 3. a) The ratio £ of effective interaction duration r and 
orbital period t or f, as a function of radius, calculated for the stan- 
dard globular cluster and normalized to one at the half light ra- 
dius Tfr. b) Dependence of the correction factor r\ on £. 



A very fast encounter of a massive black hole with a Galac- 
tic globular cluster can therefore effectively alter both the 
center of mass energy as well as the internal energy (i.e. the 
velocity dispersion: Ei n t oc a 2 ) of the cluster, whereas a very 
slow encounter only affects the center of mass velocity of the 
system. 

The orbital periods of stars in globular clusters differ by 
many orders of magnitude between stars in the core region 
and stars in the outer envelope. Thus for certain impact 
parameters b and encounter velocities dq the impulsive ap- 
proximation will be valid in the outer parts, whereas in the 
inner regions the adiabatic limit is reached. A realistic appli- 
cation of the impulsive approximation to these objects has 
to take this into account. Spitzer (1958) defined the relative 
duration of the interaction £ as being the fraction of the ab- 
solute encounter time r relative to the orbital period t or \, of 
an individual cluster star: £ = r/torb- Under the assumption 
of a harmonic cluster potential he was able to derive an ana- 
lytical expression for a correction factor rj(Q to parametrize 
the effectivity of the interaction (Av/v — > r/Av/v), 



77 = -~(L X + L y + L z 



with 



(7) 

(8.a) 
(8.b) 



MO = [C 2 ^i(C) + CMO] + [c 2 #i(0] , 
MO = [C 2 ^o(0] 2 + [C 2 ^i(0] 2 , 
MO = [CMO] 2 • (8x) 

Here the K n (Q are the modified Bessel functions of 
second rank and n-th order, the Macdonald func- 
tions. These follow for large £ the asymptotic formula 



Kn(0 = (I^C 1 ) 



(Bronstein & Se- 



mendjajew 1987). So, for small encounter times r — > (and 
thus £ — > 0) the (classical) impulsive approximation is recov- 



ered and for large durations (£ — > 00) the adiabatical limit 
is reached, which means the effectivity of the interaction 
exponentially approaches zero. 

Within the framework of this paper it is best suited to de- 
fine the encounter time r as the time interval within which 
(under idealized conditions, i.e. fixed particle location and 
straight moving perturber) 95 % of the total momentum 
transfer takes place: r = At | g5% ~ 3b/vo. The stars of 
a globular cluster move - under the assumption of circular 
orbits - with velocities v 2 = rd$/dr = GM(r)/r in the 
potential <J>(r) where M(r) = 4-7T p(r')r' 2 dr' is the mass 
within the radius r. The orbital period is then 



' ' = ~vT = V GM{r) 



(9) 



and the relative duration of the encounter follows as 
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Figure g.a shows the values of £ for the standard globular 
globular cluster (model #1, see Sect. ^) normalized to one at 
the half light radius r^. The correction factor 77 (equation ^ 
scales as plotted in Fig. ^.b. The slight increase of the inter- 
action efficiency for £ < 0.8 results from resonances between 
orbit and perturber potential. For a large interaction dura- 
tion £ 3> 1 the adiabatical limit is reached, r\ approached 
zero exponentially. 



3.3 Application to Globular Clusters — Energy 
Transfer 



To calculate the total energy input by an encounter with 
a massive black hole into a globular cluster as an extended 
object consisting of many ( 10 5 ) stars, one best applies the 
modified impulsive approximation as derived above to every 
single star and then sums up. The increase of kinetic energy 
immediately after the encounter is 



1 A 2 

Av, + -rriiAvi 



(11) 



where AK* is the change of kinetic energy of star i and m;, 
Vi and Avi are its mass, its velocity within the cluster prior 
to the encounter and the resulting velocity change, respec- 
tively. The total change of kinetic energy AK is obtained by 
the statistical average 



AK = Mv ■ Av + -A/Av 2 



(12) 



AK contains the center of mass acceleration of the whole 
cluster as well as its increase of internal energy. The later is 
the parameter of interest in the present context. It is to be 
compared with the binding energy of the system Eb ■ 

Since in this approach the stars within the cluster are con- 
sidered as fixed in space during the encounter, the potential 
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Figure 4. Relative increase of the velocity dispersion (o"/i na ; — ctinit) / <?init (solid line) and the center of mass velocity (dotted line; 
values in kms ) due to the encounter as function of impact parameter b (normalized to the tidal radius rt) and impact velocity vq of 
the black hole, computed using the modified impulsive approximation. The shaded circles denote the parameter pairs, for which N-body 
simulations were done. The total mass of the globular cluster is 1.03 X 10 Mq (model #1: standard cluster with rt = 33.71 pc) and the 
mass of the black hole is 10 6 Mq. 



energy is not altered. Due to the increase of kinetic energy, 
the globular cluster has lost its virial equilibrium and be- 
gins to relax towards a new equilibrium state. Throughout 
this process kinetic energy is converted into potential energy 
and the cluster as a whole expands. and vice versa. Thereby 
it eventually evaporates stars through its tidal barrier and in 
strong interactions (A.K > Eb) dissolves completely (Cher- 
noff & Weinberg 1990). 



3.4 Comparision with Numerical N-Body 
Simulations 

To test the accuracy of the above approximation, we have 
calculated the encounters of the standard globular clus- 



ter (model #1, see Sect, g) with a black hole of mass 
Mbh = 10 6 Mq for 14 different pairs of impact parame- 
ter b and relative velocities no using a standard Tree Code 
scheme (e.g. Barnes and Hut 1986, Hernquist 1987, 1990). 
We compare its result with the parameter study undertaken 
with the modified impulsive approximation described above. 
Figure ^ shows the relative increase of the internal velocity 
dispersion (a final — Cinit) /o~init and of the center-of-mass ve- 
locity Avcm of the cluster in kms -1 due to the encounter, 
as predicted by the impulsive approximation. The encounter 
speed varies from 10 kms -1 to 500 kms -1 , the impact pa- 
rameter - scaled in units of the tidal radius of the globular 
cluster - from (black holes penetrates through the center 
of the cluster) via 1 (the black hole just touches the cluster) 
to 10 (distant encounter). The adiabatical correction factor 
has a strong influence for velocities smaller than 250 kms -1 . 
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The amplification of the encounter effectiveness rj is due to 
resonances between perturber and the orbits of cluster stars 
(see equation (Q) or Fig. ^) and pushes the lines of equal in- 
crease of velocity dispersion to larger impact parameters b. 
In the upper half of the figure, i.e. for large encounter speeds, 
the adiabatical correction is negligible and the original hy- 



perbolic shape of the lines Aa/cr = const, is recovered. The 
center of mass motion is not affected by internal adiabatic 
effects. Therefore the lines of constant Avcm show the ex- 
pected hyperbolic behavior throughout the whole considered 
parameter range. The shaded circles denote the parameter 
pairs for which additional N-body simulations were made. 
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The results of these simulations were compared with the an- 
alytical predictions in Fig. ^. The numerically derived values 
for (a fi na i — ainit) / <Jinit and Avcm agree remarkably well 
with the analytically derived ones. Therefore the modified 
impulsive approximation is an adequate tool to determine 
the increase of the internal energy of a globular cluster due 
to such encounters and thus the probability for the cluster 
to get disrupted. 



Table 2. Average Number N of encounters between the standard 
cluster (model #1) on circular orbits at R = 5kpc, R = lOkpc 
and R = 15 kpc and black holes with masses Mbh = 10 6 Mq 
during t = 10 10 years depending on the impact parameter b = 
b/rt, normalized to the tidal radius rt of the cluster (r^ denoting 
the half mass radius). 



6 7V(5kpc) JV(lOkpc) JV(15kpc) 



4 DISRUPTION TIME SCALES FOR THE 
STANDARD CLUSTER 

We adopt in our calculations a mass model for the Milky 
Way, which contains a bulge and disk component and a dark 
isothermal halo to account for the flat rotation curve with 
Vrot — 220kms _1 . The bulge is analytically represented by a 
Plummer spheroid with total mass M B = 1.5 x 1O 1O M and 
core radius tb = 1.5 kpc. The disk component is calculated 
as double exponential with radial scale length Rd = 3.5 kpc 
and vertical scale height zn = 0.25 kpc. Its density distri- 
bution is normalized to pQ = 0.072 Mqpc -3 in the solar 
neighborhood (Allen 1973, §118). The dark halo results then 
in 



p(r) = 0.134 

pc J 



1 + 



R u 



(13) 



where Ru = 2.5 kpc denotes the core radius. In the solar 
vicinity, i.e. for r — Rq = 8.5 kpc, one obtains a dark matter 
density of p(Ro) = 0.01 Mq pc -3 . Following the considera- 
tions of Lacey & Ostriker (1985) that this dark halo consists 
entirely of massive black holes with Mbh = 10 6 Mq , there 
are about 10 of these objects in a volume of 1 kpc 3 around 
the sun. In the center of the Galaxy we would expect approx- 
imately 134 black holes per kpc 3 . The 1/r 2 - density fall-off 
resembles that of an isothermal sphere and thus one can 
derive a velocity dispersion for these objects of 



a bh = yj2ixGr 2 p(r) ~ 120km s" 



(14) 



The velocity distribution of the black holes in one direction is 
then assumed to be Gaussian with zero mean and a standard 
deviation of 120 km s -1 . 

We assume, that the globular clusters move on circular or- 
bits with v — v c ~ 220 kms~ through this background 
medium of black holes. This is a somehow simplistic ansatz, 
but adequate for the so called thick disk population of globu- 
lar clusters and a reasonable order of magnitude approxima- 
tion for the old halo clusters. Then the number of encounters 
between these black holes and an individual globular cluster 
at a galactic radius R with impact parameters in the interval 
[b, b + db] and impact velocities [v,v + dv] within a period t 
is 



N{R;b,v,t)dbdv 
= 2Mb- 
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The expression in brackets accounts for the Gaussian dis- 
tribution of black hole velocities. Thus in the Milky Way 
the total number of encounters with b < b within the time 
interval t is 



N(R;b< b,t) = 
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(16) 



The resulting numbers of encounters between black holes 
with A/bh = 10 6 Mq and a standard cluster (model #1, see 
Tab. ^) at galactocentric distances R — 5 kpc, R — 10 kpc 
and R — 15 kpc within a period of t = 10 10 years are given 
in Tab. |[ 

During its lifetime a globular cluster would experience very 
many encounters, even in the outer parts of the Galaxy. 
Clearly this bombardment would have devastating effects 
on the population of Galactic globulars, as will be quantified 
below. 

To compute the gain of kinetic energy due to interactions 
with black holes, one has to convolve the number of encoun- 
ters N(R;b,v,t) with the increase of the internal velocity 
dispersion per encounter [(a fina i - ainit) / 'ainit] (b, v) of the 
globular cluster: 
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Figure 6. Function b'Acr/cr (dotted line) and the resulting in- 
tegral (solid line) depending on the relative impact parameter 
b' = b/rt at fixed v = v c for the standard model of a globular 
cluster in the Milky Way. 



next section, which describes the Monte Carlo simulation of 
individual cluster histories within the Galaxy. 



5 MONTE CARLO SIMULATIONS OF 

INDIVIDUAL ENCOUNTER HISTORIES 

The mean free path of particles with velocities (v) in a sea 
of background particles with number density n and cross 
section E is (I) = (n'E)~ 1 — (v)t. Thus the average time r 
between two encounters with given impact parameter in the 
interval [6, b + Ab] and relative velocity v is 



1 

n'Ev 



(19) 



Using the values of the standard cluster (model #1 with 
[Aa/a](b,v) according to Fig. || and Mbh = 10 6 M©), this 
results in 



V final ® init 
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How well does the integral in equation [n] converge and 
which impact parameter dominate the result? These ques- 
tions are addressed in Fig. ^. It shows the integrand 
b' [Aa /a] (dotted line) at fixed v = v c and the resulting 
integral (solid line), as function of b' — b/rt, the impact 
parameter normalized to the cluster radius. Saturation and 
thus convergence of the integral is reached already within 
the cluster itself. This indicates that only a relatively small 
number of central encounters with b' < 1 dominates the 
total energy gain. 

The above calculation shows a large increase of the veloc- 
ity dispersion of the standard cluster over the lifetime of 
the Milky Way, e.g. by a factor of 4.7 at a galactocentri- 
cal distance R — 10 kpc. No cluster would survive such a 
bombardment. But what is the scatter around the mean 
of this value? Equation (^) was obtained by integrating 
over a smoothed background distribution of dark halo ob- 
jects. Under the assumption of black holes being the main 
constituent of the dark halo, this smoothness is no longer 
true: the distribution of matter in the halo is expected to 
be extremely coarse grained. Furthermore, the evolution of 
the globular cluster is mainly affected by encounters within 
the innermost region and these are relatively rare (as shown 
above). Which means the collision history and the result- 
ing dispersion gain is strongly determined by only a handful 
events and thus one expects huge deviations from the mean 
for different encounter histories. Therefore, even if the av- 
erage energy gain would rule out the survival of a certain 
type of globular clusters within reasonable distances from 
the galactic center, individual clusters might still be able to 
resist this destruction mechanism due to the expected large 
scatter in Aa ja as a result of different encounter histories. 
A realistic study of the survival rate within a population of 
specific globular clusters must take this effect into account 
and apply statistical methods. This shall be done in the 



The number density of halo black holes is n(R) = 
p(R)/Mbh, adopting a halo density distribution p, according 
to equation ([l3]). The interaction cross section is E = 2-nbAb 
for simple geometric arguments. Since the distribution of 
dark halo objects is assumed to be isothermal with veloc- 
ity dispersion ctbh — 120 kms -1 (see equation Q), the dis- 
tribution of encounter velocities between the halo objects 
and a globular cluster is a Gaussian with standard devia- 
tion obh around the center of mass velocity of the globular 
cluster. Assuming they are moving on circular orbits with 
v c ~ 220 kms , one has to weight the velocity entering 
equation @ with (27ro- BH )~ 1/2 exp(-|(i>--y c ) 2 /o-| H ). This 
way we obtained the typical time interval r(b, v) between 
two encounters for the whole parameter region described in 
Sect. ~~ 



3.4 



Transport theory (consult e.g. Landau & Lifschitz 1983, Vol. 
X) specifies the probability for a particle not to undergo an 
encounter in the time interval At as P(At) = e~ At ^ T . There- 
fore, the probability that a globular cluster actually does en- 
counter a black hole within At and with interaction param- 
eters (b,v), results as P(At) = 1 - P(At) = 1 - e - Ai MM)^ 
The time interval At corresponding to a given probability P 
is then At(P) = — r ln(l — P). Since encounters with differ- 
ent impact parameters b and v are statistically independent, 
one can handle each grid point (b, v) in parameter space sep- 
arately. For each pair (b, v) we throw a dice to obtain a ran- 
dom probability P and translate this random number into a 
time interval At(P) as described above. If At(P) is smaller 
than the studied total time interval to (here initially 10 10 
years), then we accept the encounter event, i.e. add Aa(b, v) 
to the velocity dispersion of the cluster and substract At(P) 
from to, leading to a remaining time t = to — At(P). This 
procedure is repeated, until the time interval At(P) to the 
next encounter becomes larger than the remaining time t. 

Figure [i] visualizes the distribution of Aa /a after 10 10 years 
for the standard cluster (model #1) in a galactic dark halo 
of 10 6 Mq black holes at three different galactic radii. The 
average dispersion gain Aa /a at these radii agrees well with 
the value obtained in the last section, using a smooth back- 
ground distribution (equation (^). However, the histogram 
shows the expected huge scatter due to the discreteness of 
the encounter events. Adopting the binding energy of the 
cluster as threshold for its disruption, we determine the frac- 
tion r/i of the clusters that get destroyed in the simulation: 
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Figure 7. Histogram of the dispersion gain Aa/a of the standard globular cluster (model #1) for a black hole mass of 10 Mq at 
R = 5, 10 and 15 kpc (10 000 encounter histories per plot). Histories that lead to Acr/cr < 1 are counted as cluste r su rvivals, rji denotes 
the fraction of these events. The average gain (Aa/a) agrees well with the analytical value obtained in equation (|18[) at these radii. 



Whenever the total kinetic energy of the cluster stars ex- 
ceeds the binding energy of the system, i.e. when Aa/a > 1, 
we assume the cluster has dissolved over the studied time 
interval. This threshold of Aa/a > 1 is a somehow crude 
approximation. We therefore have performed N-body calcu- 
lations of encounters with Aa/a^t 1, which indicate that the 
mass loss during such an event is so severe, that the cluster 
either dissolves instantaneously or the remaining (bound) 
stars relax back into an equilibrium state that is much less 
viable and will get disrupted within a few 10 8 years. Clusters 
that suffered a central encounter with Aa/a > 1.5 were al- 
ways destroyed. An event with Aa/a = 1.25 typically causes 
the cluster to loose 25 % of its stars. The energy input in a 
central encounter is not distributed evenly onto the cluster 
stars, but affects those most that are located along the tra- 
jectory of the black hole through the cluster. 75% of the stars 
remain bound and form a new globular cluster like configu- 
ration within a few relaxation time scales r rc i. The binding 
energy of this "new" cluster however is much smaller than 
that of the original one by a factor of five. Therefore the 
next encounter of similar strength (or several more distant 
encounters with the same cumulative energy input) will eas- 
ily disrupt the cluster. After all, it is save to take the value 
Aa/a = 1 as an adequate upper limit to cluster survival. 

The simulations also indicate that a series of distant (and 
therefore weak) encounters integrated over a certain period 
At have a more serious effect on a globular cluster than one 
single central encounter at the beginning of this time inter- 
val At leading to the same net energy input. In the first 
case, the excess energy is more or less evenly distributed 
onto the cluster stars thus leading to a high evaporation rate 
throughout the whole studied time interval. Whereas in the 
case of a central collision, some stars (alongside the trajec- 
tory of the black hole through the cluster) gain very high 
velocities and immediately leave the cluster taking away a 
large fraction of the energy input. The remaining stars stay 
kinematically colder and can - when the impact was not 
completely devastating in the first place - relax back into 
equilibrium, then again exhibiting the small mass loss rate 
of unperturbed dynamical evolution. 

The exemplary case, standard cluster at R = 10 kpc with 
black holes of Mbh = 10 6 Mq, reveals a extremely low 



rate of survivors within the inner 10 kpc of the Milky Way 
(Fig. ^). This analysis is extended to all models and to dif- 
ferent black hole masses in the next section. 



6 ENCOUNTER HISTORIES FOR 

DIFFERENT GLOBULAR CLUSTER 
MODELS IN THE MILKY WAY 

Determinations of the total mass in the luminous halo of 
our Galaxy lead to a value around 1.5 x 10 9 M© (Car- 
ney et. al. 1991; or e.g. M H = 9 x 10 s M in the range 4 - 
25 kpc in Suntzeff et. al. 1990). With a total mass of around 
2 x 10 7 M Q in (halo) globular clusters out to R = 40 kpc 
(Carney et. al. 1991), the cluster to halo mass ratio is 0.015 
- 0.02. The maximum allowed cluster disruption rate is ob- 
tained, assuming that all halo field stars are the debris of dis- 
solved Galactic globular clusters, implying that the Galac- 
tic globular cluster system was about a factor of 100 more 
massive when it formed more then 10 10 years ago than it 
is today. Neglecting the evolutionary differences of globular 
clusters of different masses and concentrations, the above ra- 
tio sets a very conservative lower limit to the survival rate of 
these clusters, which further can be used to obtain an upper 
limit on the mass of dark halo black holes: If the assumption 
of a certain black hole mass leads to a survival rate below 
771 = 0.01, then this mass has to be discarded. With the 
method described in the last section we calculate for each of 
the ten globular cluster models the fraction of survivors for 
black hole masses in the range 10 3 Mq to 10 7 Mq at three 
galactocentric radii, R — 5 kpc, 10 kpc and 15 kpc. The re- 
sult of these computations is plotted in Fig.'s ^, The 
dashed line indicates the acceptance threshold at rji = 0.01. 

These figures exhibit the following trends in our model: 

(i) Closer to the Galactic center the number density of 
dark halo objects increases and so does the number of en- 
counters. Thus the range where the number of surviving 
globular clusters drops to zero shifts to lower black hole 
masses. 

(ii) Within one sequence of equal concentration, the more 
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massive clusters are more viable. This is understandable, 
because the bigger the cluster is, the tighter it is bound 
at fixed concentration (i.e. the deeper is its own potential 
well) and the better can it resist the steady bombardment 
of background particles. 

(hi) Since M oc po and also M oc a 2 (according to equa- 
tion hi), the stability of the cluster increases with growing 
central density or dispersion. 

(iv) Quite unexpectedly, we find that the most stringent 
constraints come from the most centrally concentrated clus- 
ters #9 and #0 to #3. The larger c, the more extended is 
the envelope of less strongly bound stars outside the core 
and thus the geometrical cross section of the cluster. With 
an extended outer region a cluster can more effectively loose 
stars than a less centrally concentrated and thus smaller one. 
Variations of the concentration index c have the largest effect 
on the survival rate of the model clusters. It even outweighs 
the influence of the cluster mass, as can be seen in sequence 
#4, #8 and #9: The threshold of 1% survivors is reached 
for cluster #9 first, even though it has the highest mass in 
that sequence. 

The above results have to be compared with the observa- 
tion: The upper limits for black hole masses at R — 5kpc 
for the ten model clusters are the strictest constraints. 
Since our model clusters all have central densities around 
5 x 10 3 M0PC -3 and the central density is directly propor- 
tional to the total mass and thus to the depth of the poten- 
tial well of the cluster it influences the binding energy of the 
cluster and its stability. One therefore has to compare the 
model clusters with Galactic globular clusters in the same 



density range (po = 10 3 - 10 4 Mqpc -3 ), in order to obtain 
meaningful results. This is depicted in Fig. [ll]. Analogous to 
Fig.^, it shows our ten model clusters (shaded circles) and 
the Mgc _ c distribution of the Galactic globulars in the 
appropriate density range. The four thick gray lines confine 
the region of too low cluster survival rate. Globular clus- 
ters with masses Mgc and concentrations c to the left of 
these lines have a probability of less than 1% (771 < 0.01) 
to survive for 10 10 years in a Galactic dark halo consist- 
ing of black holes with masses Mbh = 5 x 10 4 M Q , 10 s M , 
2.5 x 10 5 M Q and 4 x 10 s M Q , respectively. These values are 
obtained from the calculations at R — 5kpc, since these give 
the strongest constraints for the maximum allowed black 
hole masses. Hence we have to compare them with the glob- 
ular clusters at i?^5kpc. One sees that all these clusters fall 
beyond the line M BH = 5 x 10 4 M Q . Therefore 5 x 10 4 M 
is an upper limit to the mass of compact dark halo objects. 
This limit is more than an order of magnitude lower than the 
value M BH ^ 10 6 M Q derived by Lacey & Ostriker (1985) 
in order to explain the disk heating process, that possibly 
led to the Thick Disk component in our Galaxy. Thick Disk 
formation by heating due to black holes hence can be ruled 
out. 

In our examination we concentrate on one single effect of 
Galactic environment on the evolution of globular clusters, 
the influence of massive compact objects building the dark 
halo of our Milky Way. It is beyond the scope of this pa- 
per to include further aspects of the dynamical evolution of 
Galactic globular clusters, such as disk shocking, influence 
of the Galactic bulge or loss of angular momentum of glob- 
ular cluster orbits due to dynamical friction, driving them 
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Figure 11. Distribution of mass Mqc and concentration c for the ten model clusters and for Galactic globular clusters with central 
densities in the range 10 3 Mgpc -3 < po < 10 4 Mgpc - 3 (Webbink 1985). The probability for globular clusters to survive is less then 1%, 
i.e. rji < 0.01, in the region left and above the gray lines for dark matter halos consisting of massive black holes of A4bh = 5 10 4 Mq, 
10 5 Mq, 2.5 X 10 5 Mq and 4 X 10 5 Mq, respectively. These lines were computed at R = 5kpc and have to be compared with the observed 
globulars with the same galactocentric distance (squares). 



closer to the Galactic center. Also the role of binaries and 
mass segregation in the internal evolution of star clusters are 
not considered. Whereas the latter two are negligible for our 
study, the first two can be of importance especially in the 
inner regions of the Milky Way. Neglecting these corrosive 
effects when comparing our model with observed properties 
of the Galactic population, leads to a somewhat too large 
upper mass limit for the black holes. Thus, a more realistic 
constraint on compact dark halo objects may lie below the 
values derived here. Concerning globular cluster trajectories, 
we assume for simplicity that they are moving on circular 
orbits. To examine the influence of this ansatz, we determine 
the number of encounters with compact background objects 
for globular clusters on orbits with different ellipticities in 
the Galactic potential, using the dark matter distribution 
(equation |l3| ) derived in Sect. ^| The number N of encoun- 
ters within a certain time interval t is TV = E ■ n ■ vt, with 
E being the collision cross section and n being the num- 
ber density of background objects. Calculating TV for orbits 
with the same average distance R to the Galactic center, i.e. 
the time average of R is kept fixed, but ellipticities ranging 
from zero to one, from circular to radial orbits, we see only a 
moderate increase of the encounter rate by a factor of two. 
This effect gets stronger in a more centrally concentrated 
dark matter distribution than adopted in equation (jH|) and 
vanishes in a homogeneous distribution. We conclude, that 



using circular orbits only slightly underestimates the colli- 
sion rate and thus only slightly increases the mass limit for 
black holes in the Galactic dark halo. 

Let us now relax the stringent threshold of T)i =0.01, adopt- 
ing for example a value of r)i — 0.5, i.e. if we require 50% 
of the initial clusters to survive. We then get a maximum 
permitted black hole mass which is by a factor of about 
5 smaller than before, thus pushing that value down to 
TO 4 M . Since it is very unlikely, that the halo field stars 
stem all from disrupted globular clusters and one can envi- 
sion many other and more reasonable formation scenarios, 
like an intrinsic population of halo stars or debris from ac- 
creted satellite galaxies, the maximum number of disrupted 
clusters is more likely in the range of a few tenth of per 
cent. Then the method applied in this paper leads to an 
upper limit on A/bh of 10 4 Mq . 

This number has to be compared with limits on A/bh derived 
by other authors. Moore (1993) obtained values which are 
lower by about one magnitude. He quasi-analytically com- 
puted the disruption time scale as a function of black hole 
mass for nine halo globular clusters in our Milky Way. De- 
scribing the cluster by a smooth density profile and thus ne- 
glecting the initial velocity dispersion of the cluster stars, he 
calculated the energy input per encounter. He then summed 
up encounters with randomly distributed impact parame- 
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ters and velocities until the total energy input is equal to 
the binding energy. The required time interval defines his 
disruption time td. By adjusting the black hole mass, he 
fixed td for each particular cluster to be roughly half the 
typical age of a globular cluster in the Galaxy or about half 
the Hubble time: td = 7 x 10 9 years. Thus he used only 
one encounter history to obtain his upper limit on Mbh for 
each cluster. Since each encounter history is dominated by 
only a few central impacts the scatter between different en- 
counter histories is huge (see Fig. [?]) and Moores numbers 
are plagued by error bars of at least one order of magnitude. 
We therefore conclude, that his value of Mbh ^ 1.1 x 10 3 M Q 
is still questionable. 

A different approach was taken by Rix & Lake (1993). 
They followed the original considerations of Lacey & Os- 
triker (1985) about the effects of massive dark halo objects 
on galactic disks, applying them to two nearby dwarf galax- 
ies dominated by an extended dark halo: DDO 154 and 
GR 8. DDO 154 (similar to DDO 170, Lake et al. 1990) 
results in M BH < 7 x 10 5 M Q and GR 8 in M BH 6 x 10 3 M Q . 



7 CONCLUSION 

The fact that globular cluster with masses in the range 
10 4 Mq < Mgc < 10 6 Mq and concentrations as high as 
c ~ 2 have survived for the past 10 10 years provides strong 
constraints on the masses of black holes as possible candi- 
dates for the dark matter in our Galaxy. As globular clusters 
represent only a small subpopulation of the Galactic halo, a 
large fraction of them could indeed have been disrupted by 
encounters with massive black holes forming the field star 
population of the halo, with only a few lucky candidates 
being left behind. We have calculated the survival probabil- 
ities for different model clusters in the Milky Way assuming 
a dark halo consisting entirely of massive black holes and 
demonstrated that these probabilities strongly vary for in- 
dividual encounter histories. In contrast to former models, 
we therefore propose that detailed Monte-Carlo simulations 
are required in order to determine the maximum allowed 
black hole masses. Our calculations lead to an upper mass 
limit of Mbh $5 x 10 4 Mq , which is somewhat larger than 
former estimates but on the other hand small compared to 
the black hole masses required to explain the formation of 
the Thick Disk as the result of dynamical heating by these 
massive objects. 

The calculations cannot rule out massive black holes com- 
pletely. Indeed, black holes with masses of M BH ^5xlO 4 M 
could in principle explain the observed upper concentration 
limit (c ~ 2) of globular clusters as more concentrated ob- 
jects would be disrupted. This possibility however is still 
unlikely given the fact that observations in dwarf galaxies 
indicate even smaller upper mass limits. 
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